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ABSTRACT 

Context. Abundances of iron-peak and a-elements are poorly known in Orion, and the available measurements yield 
contradictory results. 

Aims. We measure accurate and homogeneous elemental abundances of the Orion subgroups ONC and OBlb, and 
search for abundance differences across the Orion complex. 

Methods. We present flames/uves spectroscopic observations of 20 members of the ONC and OBlb. We measured 
radial velocity, veiling, effective temperature using two spectroscopic methods, and determined the chemical abundances 
of Fe, Na, Al, Si, Ca, Ti, and Ni using the code MOOG. We also performed a new consistent analysis of spectra previously 
analyzed by our group. 

Results. We find three new binaries in the ONC, two in OBlb, and three non-members in OBlb (two of them most 
likely being OBla/25 Ori members). Veiling only affects one target in the ONC, and the effective temperatures derived 
using two spectroscopic techniques agree within the errors. The ONC and OBlb are characterized by a small scatter 
in iron abundance, with mean [Fe/H] values of —0.11 ± 0.08 and —0.05 ± 0.05, respectively. We find a small scatter in 
all the other elemental abundances. We confirm that P1455 is a metal-rich star in the ONC. 

Conclusions. We conclude that the Orion metallicity is not above the solar value. The OBlb group might be slightly 
more metal-rich than the ONC; on the other hand, the two subgroups have similar almost solar abundances of iron-peak 
and a-elements with a high degree of homogeneity. 

Key words. Open clusters and associations: individual: Orion Complex - Stars: abundances - Stars: low-mass - Stars: 
pre-main sequence - Stars: late-type - Techniques: spectroscopic 



1. Introduction 

Measurements of elemental abundances in low-mass mem- 
bers of young clusters and associations represent important 
tools for addressing different issues in the field of star and 
planet formation. 

On the one hand, elemental abundance determinations 
in young associations, as for old populations, allow one to 
use chemical tagging to investigate formation scenarios and 
possible common origin of different groups. In particular, 
accurate measurements of abundances and abundance ra- 
tios in members of different regions and subgroups belong- 
ing to the same star-forming complex can unveil group- 
to-group differences and chemical enrichment. In turn, this 
would represent a signature of se quential star forma tion and 
supernova nucleosynthesis (e.g-. ICunha et al.lll998l and ref- 
erences therein). 

On the other hand, giant gas planets are preferentially 
found around old solar-type star s more metal-rich than 
the Sun (e.g.. I Johnson et al1l2010L and references therein). 
Since planets are assumed to form from circumstellar (or 
proto-planetary) disks during the pre-main sequence (PMS) 
phase, the obvious question arises about the metallicity of 
young solar analogs and what fraction of them (if any) is 
metal-rich. 
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Orion is one of the most well-studied star-forming com- 
plexes and represents an ideal laboratory for investigating 
all the stages related to the birth of stars and planetary sys- 
tems. In particular, the different ages of the four subgroups 
(~8-12 for la, ~3-6 for lb, -2-6 for lc, and < 1-3 for Id; 
Bally 2008) belonging to the OBI association, along with 
the content of dust and gas, appear to support th e idea 
of sequential star- formation scenario (jBlaauwl 1196411 . with 
the Id subgroup (the Orion Nebula Cl uster - ONC) being 
the yo ungest. As discussed in detail bv lCunha fe LambertJ 
(1992), since massive stars are a major site of nucleosyn- 
thesis, the gas from which the younger subgroups formed as 
a second generation may be contaminated by the enriched 
ejecta of the first generation of ma ssive stars (|Reeveslll972l 
119781 : iPreibisch fe Zinneckerll2006h . In this case, one would 
expect to detect different abundance pattern across the 
cluster subgroups, with the youngest regions - the ONC in 
particular - showing peculiar chemical enrichment in iron- 
peak and a-elements with respect to the older ones. Owing 
to the large number of supernovae expected to have oc- 
curred in Orion, this prediction can be tested by accurate 
and homogeneous abundance measurements in the different 
subgroups. 

Moreover, the high frequency of proto-planetary disks 
around low-mass members of the ONC suggests that plan- 
etary systems might be forming around a fraction of these 
stars. Determining the metallicity of the cluster is im- 
portant to investigating the connection between the early 
phases of planet formation and planets that formed some 
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Gyr ago and are currently detected around old solar-type 
stars. We mention that none of the star- forming regions 
(SFR s) for which a me t allicity is available is metal-rich 
( e.s.. ISantos et all l2008t iGonzklez-Hernandez et al.l 120081: 



ID'Orazi et alJl2009fl26Tffl . which is indeed puzzling. 

Several studies have been carried out in the past 
two decades designed to measure the abu ndances of the 
gas an d stars in Orion. As summarized by ID'Orazi et al.l 
(|2009() . however, these studies have yielded discrepant 
results in term of both the average metallicities of 
the different subgroups, the ONC in particular, and 
the presence of group-to-gr oup differences. For exam- 
ple, ICunha k Lambert! (|1994l) detected variations in oxy- 
gen and silicon abundances across Orion, which they in- 
terpreted a s the signature of self-enrichment. However, 
iSimon-DiazI (|2010f) analyzed high quality spectra of 13 B- 
type stars in Orion OBla,b,c,d and found a high degree 
of homogeneity, in contrast to the results of Cunha and 
collaborators. 

The most recent determination of [Fe/H] in Orion base d 
on late- type stars was performed bv ID'Orazi et al.l (|2009h . 
who analyzed a small sample of cool ONC members and one 
candidate member of the OBlb association. They inferred 
a solar metallicity for the ONC with a very small star-to- 
star dispersion ([Fc/H] =— 0.01 ± 0.04), along with hints of 
a slightly sub-solar metallicity for OBlb. While this result 
might provide support to th e sequential star formation sce- 
nario, ID'Orazi et al.l ([2009) emphasized that their results 
should be confirmed based on a larger sample of stars and 
the analysis of spectra more suitable for abundance mea- 
surements. 

Here, we present a new study of the elemental abun- 
dances in Orion, based on high-resolution spectra obtained 
with flames/uves on the Very Large Tel escope (VLT) 
Not o nly is our sample larger than that of ID'Orazi et al.l 
(|2009h . but the spectral range covered by our spectra in- 
cludes a significantly larger number of Fei and Fen lines, 
which enabled us to achieve a more secure determination 
of stellar parameters. Furthermore, for the coolest stars in 
the sample, the analysis was performed using GAIA models 
(jHauschfldt et al.l Il999t iBrott k Hauschildtj I2010D whic h 
are more appropriate than ATLAS models (|Kuruczl fl993). 
because of the inclusion of millions of molecular lines in the 
line list, a s explained in detail in Appendix IB"1 Finally, the 
sample of ID'Orazi et al.l (|2009l ) was reanalyzed. 

In Sect. [H we describe the sample, observations, and 
data reduction. The measurements of radial velocity, ef- 
fective temperature, veiling, and elemental abundance are 
given in Sects. [3] and |U The results, discussion, and conclu- 
sions are presented in Sects. [5j |H and [7] In Appendix [A] 
we give the line list, while in Appendix |B1 we describe the 
impact of model atmospheres on abundance measurements. 



2. Sample, observations, and data reduction 

2.1. The sample 



The ONC target stars were taken from iHillenbrandl (|l997f) : 
we selected stars with spectral types from late-G to early- 
M without evidence of strong accretion and, thus, spectral 
veiling. We avoided stars with large rotational velocities 
(vs'mi > 30 km/s) and known to be binaries. The total 
sample contains 10 stars. Similar criteria were applied to 



the OBlb group, where we s elected 10 late-K stars from 
the lBriceno et al.1 (|2005ll2007| ) samples. 

The ONC and OBlb samples are listed in Table [Q 
along with information from the literature. For the ONC, 
we indicate in Cols. 1-7 the star name, I magnitude, 
V — I color, spectral type, effective temp erature, lu- 
minos ity, a nd membership pr o babilit y from IHillenbrandl 
(119971 ) and iHillenbrand et al.l (I1998T ) and in Cols. 8-9 
the values of vsini (I Wolff et all 12004 ISicilia-Aguilar et all 
120051: ISantos et ail I2008D andsome notes. For OBlb, we 
list the star name, V magnitude, V — I color, spectral 
type, effective tempe r ature from spectral-type using the 
IKenvon k Hartmann (119951) sca le , lum inosity, and object 
class taken from Bricefio et al.l ([2005). In the last two 
columns of the table, we report our radial velocity mea- 
surements and comments on membership (see Sect. [3]). 

In Figs. [T] and [U we show the distribution in the sky of 
our targets. The ONC stars fall inside the main cluster, with 
the exception of P1455. The case of this star is discussed in 
Sect. 16.31 As for the OBlb targets, three of them fall close 
to the OBlb/OBla boundary defined bv I Warren k Hesserl 
([1977ft and is discussed in Sect. I3~T1 



2.2. Observations and data reduction 

The observations were obtained in 2009 with the 
Fiber Large A r ray M ulti-Element Spectrograph (flames; 
iPasauini et al.l l2002h attached to the Kueyen Telescope 
(UT2) at Paranal Observatory. Both the ONC and OBlb 
were observed using the fiber link to uves. We allocated 
fibers to a maximum of seven stars, leaving at least one 
fiber for the sky acquisition. We used the CD#3 cross- 
disperser covering the range 4770-6820 A at the resolution 
R = 47000. This setup allowed us to select around 60+9 
Fei+Fen lines, as well as spectral features of a- and Fe- 
peak elements (see Sect. 14.21) . 

The ONC was covered with two different pointings, each 
including seven stars with an overlap of four stars. For each 
pointings, we obtained three 45 min long exposures, result- 
ing in a total integration time of 2 h and 15 min. Three 
pointings were instead necessary to observe the OBlb stars. 
The pointings included four, three, and three stars. Each 
field was observed six, two, and four times, for a total ex- 
posure time of 4.5, 1.5, and 3 h. The log book of the obser- 
vations is given in Table [2] 

Data reduction was performe d using the flames/uves 
pipeline ( Modiglia ni et al.l I2004T) and the following proce- 
dure: subtraction of a master bias, order definition, extrac- 
tion of thorium-argon spectra, normalization of a master 
flat-field, extraction of the science frame, wavelength cali- 
bration of the science frame, and correction of the science 
frame for the normalized master flat-field. Sky subtraction 
was performed with the task SARITH in the IRAfQ echelle 
package using the fibers allocated to the sky. 

All the acquired spectra for each star were shifted in 
wavelength for the heliocentric correction and then coad- 
ded, after checking for possible radial velocity variations. 
The final signal-to-noise ratio (S/N) is in the range 40-200 
for the ONC stars and 30-80 for the fainter OBlb targets. 



1 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of the 
Universities for Research in Astronomy, inc. (AURA) under co- 
operative agreement with the National Science Foundation. 
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Table 1. Sample stars. For the ONC stars we list: object name, I magnitude, V — I color, spectral type, photometric effective temperature, luminosity, membership 
probability, v sin i, and notes. For the OBlb star we give: star name, V magnitude, V — I color, spectral type, effective temperature from spectral type, luminosity, 
and object class. In the last two columns our radial velocity values and comments are listed for both regions. 



ONC 



Star a ' c 


7 


V-I 


Sp. Type 


J cff 


log(L/L ) 


Mem 


v sin i 


Notes' 1 


Vrad 


Comment 6 




(mag) 






(K) 




(%) 


(km/s) 




(km/s) 




JW365 


11.70 


1.96 


K2-3 


4775 


0.82 


99 


67.6±5.8 


NIR, WTTS: 




SB2 


JW373 


13.48 


3.37 


K2-7 


4775 


0.99 


99 




HH, NIR, CTTS 




SB1 


JW589 


11.30 


2.12 


G8-M0 


5236 


1.23 


99 


14 


r 


27.1±0.8 


M 


JW601 


12.52 


1.25 


K2-5 


4775 


0.05 


98 


12.3±0.9 


WTTS 


22.0±0.4 


M 


JW641 


11.51 


1.74 


mid-G/early-K 


5236 


0.91 


99 


46.3±4.1 






SB2 


JW733 


13.14 


2.33 


M0.5 


3724 


0.04 


99 






27.4±0.1 


M 


JW868 


13.19 


1.96 


K3-5 


4581 


0.16 


99 


11.7±1.3 


NIR, HH 


24.7±0.5 


M 


JW907 


12.77 


1.60 


Kl-4 


4775 


0.17 


99 


32.8±2.2 


WTTS 


22.9±0.2 


M 


JW157 


10.15 


1.56 


KO-4 


4775 


1.19 


99 


7 


N 


28.0±0.1 


M 


P1455 


10.08 


0.76 


G0-1 


5902 


1.06 


97 


21 




21.2±0.2 


M 



OBlb 



Star b ' c 


V 
(mag) 


V-I 


Sp. Type 


J-cB 

(K) 


log(L/L ) 


Notes' 


(km/s) 


Comment 6 


CVS0118 


14.70 


1.26 


K5 


4350 


-0.43 


WTTS 


32.1±0.4 


M 


CVS0125 


14.53 


1.34 


K5 


4350 


-0.33 


WTTS 


27.3±1.4 


M 


CVS0128 


14.65 


1.50 


K6 


4205 


-0.17 


WTTS 




SB1 


CVS0129 


14.15 


1.23 


K6 


4205 


-0.17 


WTTS 




SB1 


CVS0159 


14.83 


1.85 


K7 


4060 


-0.14 


WTTS 


26. Oil. 3 


M 


CVS0161 


14.32 


1.40 


K6 


4205 


-0.19 


WTTS 


31.7±0.9 


M 


CVS0165 


13.73 


1.40 


K6 


4205 


0.15 


CTTS 


31.4±0.7 


M 


CVS056 


14.81 


1.53 


K7 


4060 


-0.33 


WTTS 


18.7±0.3 


PM 


CVS058 


14.79 


1.52 


K7 


4060 


-0.24 


CTTS 


18.7±1.0 


PM 


CVS065 


14.47 


1.35 


K6 


4205 


-0.33 


WTTS 


10.2±0.3 


NM 



a : JW= Uones fc Walker! (Il98gf) . P= lParenagol (|1954l ): b : CVSO=CIDA Varia bility Survey of Ori on (jBriceno et al.ll2005l ). 

c : 2MASS (|Cutri et al.H2003n data for all the ONC/OBlb targets; Spitzer (iRebull et all [20061 ) data at 3.6 < A < 8 /im for JW373 and JW157; Infrared (|Gezari et al.lli"999l ) 
data at 1.25 < A < 25 /im for JW157, JW373, JW589, and JW601 

d : HH=Herbiff-Haro objec t of host star; r=radio continuum source; N=N-band excess; NIR=photometric near-infrared excess. Notes from iFeigelson et all (|2002D and 
ISicilia-Aguilar et al l (|2005D . 

e SB=Spectroscopic binary (1: single-lined; 2: double-lined); M=Member; PM=Probable OBla/25 Ori member; NM=non-member of OBlb. 
f : WTTS: Weak-lined T-Tauri star; CTTS: Classical T-Tauri star. Notes from lBriceno et al.l (120051 ). 
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Fig. 1. Spatial distribution of our ONC targets (filled star s) and the sample of iD'Orazi et al.l (|2009ft re-analyzed in this 
work (empty stars). Dots represent the iHillenbrandl (|1997f ) sample with membership probability higher than 90%. The 
field is centered on the Trapezium cluster and covers an area of about 0.5° xO. 5°. The position of the star at the edge of 
the main cluster (P1455) is given. 



The coadded spectra of the stars we used for abundance 
measurements are shown in Figs. [3] and [4] 

3. Radial velocities and membership 

Since we did not acquire spectra of radial velocity (RV) 
standards, we first selected two stars to be used as tem- 
plates in the ONC and OBlb: JW601 (ONC) and CVS0118 
(OBlb). Both have high S/N spectra (~50-100), low v sin i 
(<15 km s^ 1 ), and radial velocity measurements from the 
literature. In particular, CVS0118 is the only star in our 
OBlb sample with a previous determination of RV. We 
then measured the RV from the first spectrum we ac- 
quired for both stars using the IRAF task RVIDLINES in- 
side the RV package. This task measures RVs from a line 
list and we used 25 lines in the spectral range 5800-6800 
A. For JW601, we obtain V rad =22.8±0.5 km s _1 , while 
for CVS0118 we find that V rad =32.1±0.4 km s" 1 . These 
values are in very good agreement with previous determi- 
nat ions of 22.0±0.4 km s " 1 an d 31.6 ±0.6 km s" 1 obtaine d 
by ISicilia-Aguilar et al.l (|2005f ) and IBriceno et all (|2007l) . 
respectively. 

Considering JW601 and CVS0118 as templates, we 
measured the heliocentric RV of all the ONC/OBlb stars 
using the task FXCOR of the IRAF package RV. To 
take advantage of the wide spectral coverage offered by 
FLAMES/UVES, we cross-correlated all the spectral range of 
our targets with the template, excluding the regions con- 
taminated by broad emission lines (e.g., Ha) or by promi- 



nent telluric features (e.g., the O2 series at A ~ 6275 A). 
To determine in the most reliable way the centroids of 
the cross-correlation function (CCF) peaks, we adopted 
Gaussian fits. The errors in the RV values were computed 
using a procedure inside the FXCOR task that considers the 
fitt ed peak height and the antisymmetric noise as described 
by iTonrv fc Davisl (|l979ft . Since we acquired several spec- 
tra per stars, we computed an average RV value for all our 
targets that did not display evidence of binarity. The aver- 
age RV values for most probably single stars are listed in 
Tabled! 

3.1. Membership 

To confirm the membership of our stars, the distributions 
of average RV measurements for the ONC are sh o wn in 
Fig.[5J along with the value derived bv lBiazzo et all (12009ft 
for ^100 very low-mass members. Seven out of 10 stars of 
ONC are confirmed as members and most probably single 
stars with a mean RV of 24.9±2.6 k m s" 1 in good agree- 
ment with previous de terminations (|Sicilia-Aguilar et al.l 
120051: iBiazzo et al.ll2009ft . The three exceptions are JW373, 
JW365, and JW641, which exhibit a double/triple-peaked 
CCF. We, thus, classify these stars as, previously unidenti- 
fied, binaries. Two of them (JW365 and JW641) also have 
a double-lined system (SB2) and are therefore discarded 
from further analysis. 

For OBab, similarl y, we compare our sa mple with the 
distribution found by IBriceno et all (|2007ft based on 30 
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Fig. 2. Spatial distribution of our OBlb stars (dots). The field covers an area of ^3° x 1°. The position of e Ori (one of the 
Orion Belt stars) is also shown. The solid line outlines the boundary between Orion OBla and OBlb (W arren &: Hesseil 
[1971 . 



Table 2. Log of the observations. 



a (J2000) 


S (J2000) 


Date 


UT 


^cxp 


# 


(h:m:s) 


(°: ' : ") 


(d/m/y) 


(h:m:s) 


00 


(stars) 


ONC 


05:35:12 


-05:23:00 


23/01/2009 


00:56:05 


2775 


7 


05:35:12 


-05:23:00 


23/01/2009 


01:44:32 


2775 


7 


05:35:12 


-05:23:00 


23/01/2009 


02:45:16 


2775 


7 


05:35:16 


-05:24:20 


28/01/2009 


02:04:19 


2775 


7 


05:35:16 


-05:24:20 


30/01/2009 


01:56:17 


2775 


7 


05:35:16 


-05:24:20 


30/01/2009 


02:44:12 


2775 


7 


OBlb 


05:29:27 


-01:31:60 


22/01/2009 


00:54:00 


2775 


3 


05:29:27 


-01:31:60 


31/01/2009 


00:59:31 


2775 


3 


05:29:27 


-01:31:60 


31/01/2009 


01:47:38 


2775 


3 


05:29:27 


-01:31:60 


01/02/2009 


00:53:31 


2775 


3 


05:33:48 


-01:25:59 


28/01/2009 


01:02:48 


2775 


4 


05:33:48 


-01:25:59 


14/02/2009 


01:41:02 


2775 


4 


05:33:48 


-01:25:59 


18/02/2009 


00:34:47 


2775 


4 


05:33:48 


-01:25:59 


18/02/2009 


01:22:50 


2775 


4 


05:33:48 


-01:25:59 


24/02/2009 


00:31:30 


2775 


4 


05:33:48 


-01:25:59 


24/02/2009 


01:32:14 


2775 


4 


05:38:35 


-01:11:60 


16/01/2009 


01:20:37 


2775 


3 


05:38:35 


-01:11:60 


21/01/2009 


02:37:49 


2775 


3 



OBlb targets from which they derived a mean RV of 
30.1 ± 1.9 km s _1 . Five sample stars are confirmed sin- 
gle members, with an average RV of 29.9±3.0 km s _1 , 



while three are non-members and two (CVS0129 and 
CVS0128) have double/triple-peaked CCFs. Among the 
non-members, CVS056 and CVS058 have R Vs close to the 
OB la /25 Ori region (mean RV ~20 km s~ 1 HBricefio et ah! 
2005), as also implied by their spatial location close to the 
OBlb/OBla boundary (Fig.©. 



4. Abundance analysis 

To summarize, abundances were obtained for all the ONC 
and OBlb single members or single-lined binaries, with the 
exception of rapid rotators (JW907, CVS0128, CVS0129, 
and CVS0161) and the very cool star JW733 for which 
measuring abundances from line equivalent widths (EW) 
is not suitable. Since CVS056 seems to belong to the 
OBla/25 Ori subgroup, we measured its iron abundance 
to gauge the properties of this region. On the other hand, 
we have not been able to derive the metallicity of CVS058 
because of its rapid rotation. 

The an alysis was per formed using the 2002 version 
of MOOG (jSnedenl Il973h that assumes local thermody- 
namic equilibrium (LTE) and where the radiative and 
Stark broadening are treated in a standard w ay. For col- 
lisional broadening, we used th elUnsoldl (ll955h approxima- 
tion. Both iKurucd (|1993l ) and iBrott fc Hauschildt] l|2010j) 
grids of plane-parallel model atmospheres were used for 
stars warmer and cooler than ~ 4400 K, respectively (see 
Appendix iBj) . This is the major ch ange introduced by o ur- 
selves with respect to the study of iD'Orazi et al.l (|2009h . 



5 



K. Biazzo et al.: Chemical abundances in the ONC and OBlb 



0) 
N 

o 



o 



1 .0 

0.8 

0.6 

0.4 

0.2 

1 .0 

0.8 

0.6 

0.4 

0.2 
1.1 
1 .0 
0.9 
0.8 
0.7 
0.6 

0.5 
0.4 

1 .0 
0.8 
0.6 



" Fe ' V ' " VI Ll t ! V I Fel VI _ 




Fel VI 

Til 

JW8B8: T„ = 4100 K_ 



Fel + Nil KM Ori: T ef) = 4700 






JW589: T ,, = 4750 K_ 





JW601: T,, = 4800 K_ 



ft 



JW373: T„„ = 5200 K_ 




P1455: T.„=5950 K_ 



20 



'avelength (A 



Fig. 3. Portion of spectra in the 6220-6270 A wavelength range of the ONC stars for which we measured the abundances. 
Different features used for temperature determination using line-depth ratios and abundance measurements are indicated. 



4.1. Spectral veiling 

We estimated the amount of veiling that affects the spec- 
tra of our stars fo llowing the procedure described by 
iD'Orazi et al.l (|2009h . In particular, we selected the nine 
lines in their list included in our spectral range. We then 
compared the equivalent widths of these lines to those mea- 
sured in the spectra of 16 members of the open clusters 
IC 2602 and IC 239 1, which are old enough (30-50 Myr; 
iRandich et all l200lh to ensure that their spectra are not 



affected by veiling. These IC stars have effective temper- 
atures similar to the ONC/OBlb targets (~4300-5800 K) 
and their spectra have a resolution close to ours. For each 
line, we checked whether any dependence of the EW on 
effective temperature was present; since for most lines we 
found a weak trend within an interval of 1000 K, we de- 
cided to bin the whole temperature range in 500 K steps 
and derive the mean EWs of the IC stars inside each bin. 
For two lines, namely Cai 6102.7 A and Cai 6122.2 A, we 
found significant trends at all temperatures, hence we used 
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Fig. 4. As in Fig. [3] but for the OBlb sample. 



the EW value of the IC cluster star with temperature closer 
to that of our ONC/OBlb star. For each line, we then ob- 



tained the veiling as nine = 



EW IC 



1. The mean 



EWoNC/OBlb 

veiling < r > was then computed as the average of all ry u 
values. 



By applying this method, we determined a veiling value 
consistent with zero for all the stars, with the exception of 
JW373 (r = 0.128±0.080) in the ONC, which is indeed clas- 
sified as CTTS (Table [lj . We thus corrected the measured 



EWs of all the lines using the relationship between the true 
EW and the measured one: £W t rue = EW mc&s (l+ < r >). 



4.2. Line list, solar analysis, and EWs 



We adopted the line list of lRandich et al. (12006ft i ntegr ated 
with lines from the list of iD'Orazi fc Randier] (|2009t ) in- 
cluded in our wavelength range. We refer to both papers 
for details on atomic parameters and their sources. The 
line list is given in Table [A~T1 



7 



K. Biazzo et al.: Chemical abundances in the ONC and OBlb 



' 1 



■ 



, This work (ONC) 
I This work (OBI b) 
. Biozzo et ol. (2009) 
. Briceno et ol. (2007) 




Fig. 5. Radial velocity distribution for the most probable 
single stars in the ONC (solid lines) and OBlb (dashed 
lines) samples. The thick solid and dashed lines refer to this 
work. The th i n solid line represents the Gaussian fit to the 
iBiazzo et al.l (|2009() sample obtained for 96 ONC targets 
with a mean V^ad =24.87 km s" 1 (ov ra d =2.74 km s -1 ). 
The thin dashed line represents the histogram obtained by 
iBriceno et all (|2007t ) for 30 OBlb stars with mean V la d — 
30.1 ±1.9 km s" 1 . 



As usually done, our analysis was performed dif- 
ferentially with resp ect to the Sun. We analyzed the 
iRandich et al.l (|2006l ) solar spectrum obtained with uves, 
using our combined line list and their solar parameters 
(T cS = 5770 K, log g = 4.44, £ = 1.1 km s 4 ). We ob- 
tained logn(Fe) = 7.52 ± 0.02 for iKurucd (ll993h models 
and logn(Fe) = 7.51 ± 0.02 for lBrott fe Hauschildtl (pOlOl) 
models. The results for all t he elements are given in Ta ble [5] 
together with those g iven bv I Anders fe Grevessd (|l989l ) and 
Aspl und et al.l (|2009l ) . We caveat that the latter values were 
obtained using 3D models. Table [3] shows good agreement 
between the two model atmospheres and between our and 
literature values. 

The EWs of the target stars were measured by means 
of a direct integration or Gaussian fitting procedure using 
the IRAF SPLOT task. Very strong lines (EW It 150 mA), 
which are most affected by the treatment of damping, were 
excluded from the list; furthermore, a 2-er clipping was ap- 
plied to the Fel list before determining stellar parameters 
and iron abundance. The abundance of the other elements 
was derived using the same criteria. 



Spectroscopic temperatures were then derived in two 
different ways. It has been demonstrated that line-depth 
ratios (LDRs) are powerful tools for measuring the effective 
temperature with a precision as small as 10-50 K for spectra 
with S/N > 100 (iGrav -fe JohansorJ ll991l iKovtvukh etaTl 
120061 IB iazzo et al.ll2007l and references therein). The pre- 
cision of this method can be improved by averaging the 
results from several line pairs. To develop appropriate T e g- 
LDR calibrations, we considered the synthetic stellar l i- 
brary described and made available bv lCoelho et al.l ([2005). 
These spectra are sampled at 0.02 A, range from the near- 
ultraviolet (300 nm) to the near-infrared (1.8 /im), and 
cover the following grid of parameters: 3500< T g <7000 
K, 0.0< log g <5.0, -2.5<[Fe/H]<+0.5, a-enhancement 
[a/Fe]=0.0, 0.4 and microturbulent velocity £=1.0, 1.8, 2.5 
km s^ 1 . 

Fo llowing the prescriptions given by IBiazzo et al.l 
(120071 ) . we used their line list in the 6190-6280 A spec- 
tral range and their line pairs, namely 15 (see their Tables 
1 and 2; we refer to that paper for a detailed explanation 
and justification of the line list and line pairs selected). 
We then measured their LDRs according to the guidelines 
of ICatalano et al.1 (|2002ft . and developed T ff-LDR calibra- 
tions after degrading the synthetic spectra to our resolu- 
tion. The calibrations were created for 3.0< logg <4.5, 
4000< T e g <6500 K, and usini=0 km s _1 , considering the 
synthetic spectra at [Fe/H]=0.0, [a/Fe]=0.0, and £=1.0 km 
s" 1 . 

In the end, the effective temperatures obtained from all 
the useful LDRs for each ONC/OBlb target were averaged 
to increase the precision of the temperature determination. 
The values (T^ s ) are listed in Table H] and plotted in Figs.© 
andE 

As commonly done, effective temperatures were also de- 
termined by imposing the condition that the Fe I abundance 
does not depend on the excitation potential of the lines. 
These temperatures are defined T c s ff (TableHJ) and represent 
the adopted values for the abundance analysis (Sect. [3]). 

We first note that the two values T^ s and T G s ff closely 
agree (see Figs. El and [7]), with the only exception of JW373, 
which is, as mentioned, a probable binary affected by veiling 
that has a rather low S/N spectrum. Moreover, in Fig.[5]thc 
agreement of both spectroscopic temperatures with T^ s is 
also good, with the exc eption of JW 8 68 an d JW589, which 
we find cooler than the iHillenbrandl (|1997ft values on aver- 
age by ~550 K. For JW733, we can give only an upper limit 
because T ff-LDR calibrations are suitable for temperatures 



Tabl e 3. Comparis on b etween solar abund a nces derived 
using |Kuruc3 (jl993ft and lBrott fe Hauschildtl (120101) model 
atmos pher es. The standard value s from Anders fe Grevessel 
(|1989f ) and lAsplund et all (|2009l ) are also listed. 



4.3. Stellar parameters 
4.3.1. Effective temperatures 

Photo metric temperatures were taken from IHillenbrandl 
(11997ft for the ONC (T c p ff in Table Q}, while for 
the OBlb sample we converted the spectral-types of 
Briceno et al.l (12005ft to tem peratures using the scale of 
Kenvon fe Hartmannl (jl995ft {T^ in Table [I}. 



Element 


log 7i ATLAS 


log TlGAIA 


log 71AG89 


log nAS09 


Nal 


6.31±0.03 


6.29±0.03 


6.33 


6.24±0.04 


All 


6.48±0.03 


6.47±0.03 


6.47 


6.45±0.03 


Si I 


7.56±0.03 


7.53±0.03 


7.55 


7.51±0.03 


Cai 


6.35±0.03 


6.34±0.03 


6.36 


6.34±0.04 


Til 


4.97±0.02 


4.97±0.02 


4.99 


4.95±0.05 


Fel 


7.52±0.02 


7.51±0.02 


7.52 


7.50±0.04 


Nil 


6.26±0.03 


6.24±0.02 


6.25 


6.22±0.04 
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Table 4. Astrophysical parameters and chemical abundances derived from our analysis. 



Star 


off 




log 3 s 


logg F 




[Fel/H] 


[Fen/H] 


[Na/Fc] 


[Al/Fe] [Si/Fc] 


[Ca/Fe] 


[Ti i/Fc] 


[Tin/Fc] 


<[Ti/Fc]> 


[Ni/Fe] 




(K) 


(K) 






(km/s) 




















ONC 


JW373 


4924±90 


5200 


3.2 


3.2 


1.8 


-0.15±0.17(18) 


-0.18±0.25(4) 


0.18±0.20 


0.13±0.23 0.10±0.20(1) 


0.04±0.18(3) 


0.00±0.17(2) 






-0.07±0.24(3) 


JW589 


4620±44 


4750 


3.5 


3.5 


1.9 


-0.16±0.01(29) 


-0.16±0.07(3) 


0.17±0.10 


0.21±0.20 0.14±0.11(3) 


0.02±0.09(2) 


0.00±0.10(5) 


0.14±0.09 


0.07±0.13 


-0.08±0.13(9) 


JW601 


4821±60 


4800 


1.0 


i.l 


2.3 


-0.14±0.10(38) 


-0.12±0.13(5) 


0.11±0.09 


0.14±0.07 0.14±0.12(5) 


-0.02±0.08(2) 


-0.08±0.08(7) 


0.18±0.07 


0.05±0.11 


-0.03±0.07(13) 


JW733 


<4159±121 






3.1 






















JW868 


4033±84 


4100 




3.8 


1.0 


-0.15±0.20(32) 




-0.14±0.23 


-0.03±0.24 0.15±0.36(2) 


-0.19±0.22(3) 


-0.22±0.21(4) 


0.40±0.21 


0.09±0.29 


0.03±0.23(7) 


JW157 


4630±33 


4700 


3.1 


3.0 


1.9 


-0.15±0.08(40) 


-0.12±0.05(7) 


0.12±0.09 


0.12±0.14 0.14±0.13(4) 


0.02±0.13(4) 


-0.05±0.11(8) 


0.20±0.10 


0.08±0.15 


-0.10±0.13(19) 


P1455 


5826±88 


5950 


4.1 


3.8 


1.8 


0.11±0.09(32) 


0.11±0.07(5) 


0.03±0.10 


0.08±0.11 0.01±0.13(2) 


0.00±0.13(3) 


0.02±0.13(6) 


-0.13±0.10 


-0.06±0.16 


-0.10±0.13(12) 


ONC a 












-0.11±0.11 


-0.09±0.12 


0.12±0.06 


0.14±0.05 0.11±0.05 


0.01±0.02 






0.05±0.06 


-0.06±0.05 


ONC 6 












-0.15±0.01 


-0.15±0.03 


0.15±0.04 


0.15±0.04 0.13±0.02 


0.02±0.03 






0.07±0.01 


-0.05±0.05 


ONC c 












-0.11±0.08 


















ONC d 












-0.13±0.03 


















OBlb 


CVSOH8 


4222±45 


4300 




4.3 


1.5 


0.01±0. 12(39) 




-0.19±0.12 


-0.10±0.12 0.11±0.13(2) 


-0.19±0.16(3) 


-0.30±0.14(6) 


0.43±0.14 


0.07±0.20 


-0.03±0.16(19) 


CVS0125 


4226±63 


4300 




4.2 


1.7 


-0.03±0.12(42) 




-0.15±0.12 


-0.04±0.12 0.15±0.16(2) 


-0.16±0.12(2) 


-0.33±0.13(6) 


0.47±0.13 


0.07±0.19 


0.03±0.15(11) 


CVS0159 


4139±116 


4000 




3.9 


1.9 


-0.12±0.23(21) 




-0.09±0.24 


-0.09±0.23 


-0.10±0.23(1) 


-0.36±0.32(3) 






-0.01±0.23(2) 


CVS0165 


4191±127 


4250 




3.7 


1.4 


-0.06±0.06(22) 




-0.18±0.07 


0.02±0.06 


-0.06±0.12(2) 


-0.22±0.21(5) 






0.04±0.31(7) 


CVS056 


4002±38 


4000 




4.0 


1.2 


-0.08±0.15(33) 


















CVS058 


4047±42 






1.0 






















CVS065 


4180±79 


4100 




4.2 


1.6 




















OBlb 












-0.05±0.05 






0.12±0.01 








0.07±0.01 


0.01±0.03 



a Mean abundances with P1455. 
b Mean abundances without P1455. 

c Mean abundances with P1455 and the D'Orazi et al. re-analysis (see Table [SJ. 
d Mean abundances without P1455 and the D'Orazi et al. re-analysis (see Table[5]). 
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higher than 4000 K. In Fig. [JJ the agreement with T^g is 
also good within the errors. 

4.3.2. Microturbulence velocities and surface gravities 

The microturbulence velocity £ was determined by impos- 
ing that the Fel abundance is independent on the line 
equivalent widths. The initial microturbulence velocity was 
set to 1.5 kms -1 . The values of £ are listed in Table |U We 
n ote that our det e rmina tions are typically higher than those 
of iD'Orazi et all (|2009f ). We comment on this in Sect.[ 
At variance with lD'Orazi et al.1 (|200£ 



who did not have 
enough Fc II features in their spectral range, we were able to 
estimate the surface gravity logg by imposing the Fei/Feii 
ionization equilibrium (log g in Table [4]) for all the ONC 
stars with the exception of JW868, where no Fell line was 
found. For the stars in the ONC, the initial log g was ob- 
tained from the relation between mass, luminosity, and tem- 
perature (log g = AAA + log M + 41og(T off /5770) - logL, 
labeled as log g p in Table H]) taking as astroph ysical param- 
eters the values given bv lHillenbrandl (|l997h . We verified 
that t he effect on the g ravity of considering T G s ff instead 



of the iHillenbrandl (|1997 ) temperature is negligible for our 
accuracy. We note that the agreement between logg s and 
\ogg p is quite good, the difference being at most 0.30 dex. 

On the other hand, since the stars in OBlb are cooler 
than those in the ONC, we did not find any Fe n lines, with 
the only exception of a couple of lines in CVS0118 and 
CVS0125. Thus, we decided to fix the surface gravity to 
the values obtained using the relation between M, L, and 
T R ff, where the astr ophysical parameters were taken from 
iBriceno et all (|2005t) . 

Fin ally, we remeasured [Fe/H] for the ID'Orazi et al.l 
(2009) stars applying the same method used here. We mea- 
sured the line EWs and derived the atmospheric parame- 
ters and iron abundances using ATLAS and GAIA models 
for stars with effective temperatures higher and lower than 
4400 K, respectively. 



4.4. Errors 

Derived abundances are affected by random (internal) and 
systematic (external) errors. 

Sources of internal errors include uncertainties in atomic 
and stellar parameters, measured equivalent widths, and 
the veiling determination. 

Uncertainties in atomic parameters, such as the transi- 
tion probability (log gf ), should cancel out, since our anal- 
ysis is carried out differentially with respect to the Sun. 

Errors due to uncertainties in stellar parameters (T c f{, 
£, log g) were estimated first by assessing errors in stellar 
parameters themselves and then by varying each parame- 
ter separately, while keeping the other two unchanged. We 
found that variations in T e g larger than 60 K would in- 
troduce spurious trends in log n(Fe) versus the excitation 
potential (x), while variations in £ larger than 0.2 km s^ 1 
would result in significant trends of log n(Fe) versus EW, 
and variations in logg larger than 0.2 dex would lead to 
differences between logn(Fei) and logn(Feii) larger than 
0.05 dex. The above values were thus assumed as uncer- 
tainties in stellar parameters. Errors in abundances (both 
[Fe/H] and [X/H]) due to uncertainties in stellar parame- 



6000 




4000 



4500 5000 
Te/ (K) 



5500 



6000 



Fig. 6. Spectroscopic effective temperatures versu s the 
photometric values obtained by IHillenbrandl (|1997l ). Our 
T^q and T c s ff values obtained using the LDR method and 
MOOG code are shown as filled and empty stars, respec- 
tively. The arrow represents an upper limit. T he error bars 
on th e x-axis refer to the values given by IHillenbrandl 
(|l997f ). while those on the y-axis refer to the uncertain- 
ties in Tjg only, while typical uncertainties in T c s ff are 60 K 
(see text). 



ters are summarized in Table [6] for one of the coolest and 
the warmest stars in our ONC/OBlb samples. 

As for the errors due to uncertainties in EWs, our spec- 
tra are characterized by different S/N ratios and it is 



4 600 



4400 - 



* 4200 - 



s 4000 - 



3800 - 




3800 



4000 



4200 



4400 



4600 



Te„ ST (K) 



Fig. 7. Our spectroscopic effective temperatures versus the 
values obtained by c onverting the spectral types given by 
Briceno et al. (120051) i nto te mperature using the tables of 



iKenvon k, Hartmannl (|1995l ). Our T^ s and T e s ff values ob 
tained using the LDR method and MOOG code are shown 
as filled and empty circles, respectiv ely. The error bars o n 
the cc-axis refer to the values given bv lBriceno et al.l (|2005l ). 
while those on y-axis refer to the uncertainties in T^ s only, 
while typical uncertainties in T c s ff are 60 K (see text). Note 
the different temperature range from Fig. [5] 
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Table 5. Our reanalysis of the lD'Orazi et all (|2009l ) sample (left) and lD'Orazi et all (|2009f) outputs (right). The number 
of lines employed for the iron abundance measurements are enclosed in parentheses. 



Name 




logs' 




[Fe/H] our 


T e g 


logs' 




^p^yjjjD'Orazi et al. 




(K) 




(km/s) 




(K) 




(km/s) 




268 


4300 


3.9 


1.8 


-0.16±0.13(26) 


4300 


3.9 


1.6 


-0.08±0.16(25) 


683 


4400 


3.3 


1.8 


-0.13±0.06(21) 


4250 


3.3 


1.6 


0.00±0.11(22) 


487 


4400 


3.9 


1.6 


-0.06±0.07(23) 


4300 


3.9 


1.7 


-0.07±0.08(18) 


223a 


4550 


3.8 


2.0 


-0.10±0.09(16) 


4450 


3.8 


1.6 


0.03±0.06(18) 


673 


4700 


4.0 


1.7 


-0.10±0.14(19) 


4700 


4.0 


1.5 


0.00±0.10(18) 



not possible a priori to estimate a typical error in EW. 
However, random errors in EW are well represented by 
the standard deviation around the mean abundance deter- 
mined from all the lines. These errors are listed in Tabled 
where uncertainties in [X/Fe] were obtained by quadrati- 
cally adding the [Fe/H] error and the [X/H] error. When 
only one line was measured, the error in [X/H] is the stan- 
dard deviation of three independent EW measurements. 
The number of lines employed for the abundance analy- 
sis is listed in Table @] in parentheses, with the exception of 
those elements (sodium, aluminium, and ionized titanium) 
where only one or two lines were used. 

Finally, as described in Sect. 14. 1[ the veiling of all 
the sample stars is negligible, with the only exception of 
JW373. For this star, we estimate that the veiling contri- 
bution to the abundance error is on the order of 0.05-0.11 
dex, depending on the element. 

The greatest contribution to the external or systematic 
errors originates in the abundance scale, which is mainly 
affected by the choice of the model atmospheres. This er- 
ror source is discussed in Appendix [B] Here, we emphasize 
that we used for both the Orion subgroups the same pro- 
cedure, instrument set-up, and prescriptions to derive the 
abundances. As a consequence, differences in abundance 
between ONC and OBlb should not be influenced by sys- 
tematic errors linked to the abundance scale. 



5. Results 

5.1. Metallicity 

Our final abundances are listed in Table 



3 1) The comparison 

between the new and the lD'Orazi"e t al. (2009) iron abun- 



dances is given in Tableland shown in Fig. [8] We comment 
on any differences in Sect. 15.31 In Figs. HO [HB and [TTJ we 
present [Fe/H], [X/H], and [X/Fe] as a function of T e g for 
the ONC and OBlb s amples. Figure p inc ludes both our 
targets and those from lD'Orazi et al~l (|2009f ). 

Figure [HI shows that, with the only exception of P1455 
in the ONC, we do not find a major star-to-star difference 
in [Fe/H], with a remarkable agreement between the iron 
abundance of the present sample and that of D'Orazi et 
al. Excluding the D'Orazi et al. sample, the mean ONC 
[Fe/H] is -0.15 ±0.01 (without P1455) and -0.11 ± 0.11 
(with P1455). Including those stars, we find similar val- 
ues, namely —0.13 ±0.03 and —0.11 ±0.08, respectively. As 
for OBlb, the star-to-star difference is minimal, consider- 
ing the rather large uncertainty that affects the measure- 
ment of the coolest star. The mean for OBlb is [Fe/H] = 
-0.05 ± 0.05, i.e. almost 0.1 dex above the ONC, although 
marginally consistent with it. 



0.0 



x 



-0.2 



-0.3 -0.2 -0.1 0.0 0. 
[Fe/Hf 0razi ot aL 



Fig. 8. Our iron abundance measurements of the 
ID'Orazi et all (|2009D sample versus their outputs. 



Moreover, Fig. [9] does not reveal any trend between 
[Fe/H] and effective temperature, with the exception again 
of P1455, the warmest star of the sample. Therefore, we 
believe that the difference between the [Fe/H] values of the 
ONC and OBlb is not due to systematic effects related to 
the effective temperature. 

Finally, we mention that for the likely OBla/25 Ori 
member (CVS056) we find [Fe/H]=-0.08 ± 0.15, a value 
very close to that of the ONC and OBlb (see Fig. p. 

5.1.1. [Fe/H] difference between ONC and OBlb? 

As mentioned in the Introduction, ID'Orazi et al.l {2009) 
found that the only star of OBlb ( HD294297 ) was ~ 0.1 dex 
more metal-poor than the ONC. ICaballerol (pOlOl) demon- 
strated that it is a non-member of the association. Our new 
results show that the average [Fe/H] of OBlb is slightly 
higher than that of ONC. The question then arises of 
whether the ONC is intrinsically more metal-poor than 
OBlb or there are systematic effects in the analysis. 

We note that NLTE effects have indeed been found to 
be important for cool stars with relatively high gravity , 
leading to an over estimate of the [Fe/H] ([Takedal 120081 : 
ISchuler et all I2010D . In Fig. [H we show [Fe/H] versus 
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Fig. 10. [X/H] versus T s ff . Stars and circles represent ONC and OBlb targets, respectively. Mean ONC and OBlb values 
and ±1(7 bars are shown as dashed areas, 45° and 215° oriented, respectively. For the Na, Al, Ca, and Ti abundances, we 
show the mean values obtained from all the ONC targets with T e s ff > 4500 K, while for the Si, Ni, and mean Ti abundances 
the average was computed considering all the ONC / OBlb stars. The horizontal error bar in all plots represents the typical 
uncertainty in T c s ff . 
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Fig. 11. [X/Fe] versus T e s ff . Symbols as in Fig. [TU1 
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Table 6. Internal errors in abundance determination due to uncertainties in stellar parameters for one of the coolest star 
(namely, CVS0118) and for the warmest star (namely, P1455) in our samples. Numbers refer to the differences between 
the abundances obtained with and without the uncertainties in stellar parameters. 



CVS0118 


T cH = 4300 K 


log 3 = 4.3 




£ = 1.5 km/s 


A 


AT cfi = -/ + 60 K 


Alogg = -/ + 0.2 




= -/ + 0.2 km/s 


[Fe i/HJ 


0.02/-0.02 


-0.02/0.02 




0.05/-0.05 


[Na/Fe] 


-0.06/0.07 


0.07/-0.07 




-0.03/0.03 


[Al/Fe] 


-0.03/0.04 


0.03/-0.02 




-0.03/0.03 


[Si/Fe] 


0.02/-0.06 


-0.06/0.01 




-0.07/0.02 


[Ca/Fe] 


-0.07/0.07 


0.07/-0.08 




-0.01/0.00 


[Tii/Fe] 


-0.07/0.08 


0.03/-0.03 




0.05/-0.04 


[Tili/Fe] 


0.03/-0.03 


-0.08/0.07 




-0.02/0.01 


[Ni/Fe] 


0.01/-0.01 


-0.03/0.03 




-0.02/0.02 


P1455 


Teff = 5950 K 


\ogg = 4.1 




£ = 1.8 km/s 




AT cfi = -/ + 60 K 


Alogg = -/ + 0.2 




= -/ + 0.2 km/s 


[Fei/HJ 


-0.05/0.03 


0.00/-0.02 




0.05/-0.06 


[Fen/H] 
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[Na/Fe] 
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[Al/Fe] 
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[Si/Fe] 
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[Ca/Fe] 
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[Ni/Fe] 
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log g for our ONC/OBlb stars cooler than 4500 K; the 
figure shows a moderate increase in the abundance with 
log<? > 4.0. Although we cannot quantitatively estimate 
the amount of NLTE effects, we suggest that they might ac- 
count for the difference between the ONC and OBlb metal- 
licities. 



5.2. Other elements 

In Figs. rrU] and [m we plot [X/H] and [X/Fe] as a func- 
tion of T c ff for the different elements derived in this study. 
These figures show that no trends are present for silicon and 
nickel, with the only exception of P1455, which is a 'rich' 
star in terms of all elements, included iron (see Sect. 16. 3[) . Its 
values of [X/Fe] are instead close to the solar ones fFig. fTTj) . 
In contrast, Na, Al, Ca, and Til are lower for stars cooler 
than ^4500 K both when considering [X/H] and when con- 
sidering [X/Fe]. 



★ ONC: Our sample 
ftONC: D'Orazi et al. sample 

• OBlb 
A OB1 a 
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Fig. 9. [Fe/H] versus T c s ff for t he OBlb targe t s and the 
ONC stars observed by us and bv lD'Orazi et all (|2009t ) and 
reanalyzed by ourselves. The position of CVS056 (OBla) 
is also shown. The horizontal bar represents the typical un- 
certainty in T c s ff . 



0.2 



-0.2 



-0.3 



3.0 



3.2 



3.4 



3.6 3.8 

log g 



4.2 



Fig. 12. [Fe/H] versus \ogg for the coolest stars (T e g < 
4500 K) of our and D'Orazi et al. re-analyzed samples 
(stars: ONC; circles: OBlb). 
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Similar results have be e n found by several authors (e.g., 
Schuler et al.l 120031 [2001 lYong et "all 12004: iBeirao et al.1 
20051 iGilli et al.l 120061: ID'Orazi fc Randichl l2009Tl and as- 
cribed to NLTE effects. In particular. ID'Orazi fc Randichl 
(|2009T) in their analysis of the young clusters IC 2602 
and IC 2391 show that, while the [Til/Fe] ratio decreases 
with decreasing temperature, [Tin/Fe] increases, suggest- 
ing that over-ionization is at work. Owing to their young 
age, cluster stars are characterized by enhanced levels of 
chromospheric activity and are more affected by NLTE 
over-ionization. As suggested by the aforementioned au- 
thors, deviations from LTE are possibly responsible for the 
effects seen in Fig. QTJ for Na, Al, Ca, and Ti elements. On 
the other hand, Ni and Si lines have higher ionization po- 
tentials (7.63 and 8.15 eV, respectively) than Na, Al, Ca, 
and Ti (~ 5.14 — 6.82 eV) and are thus less affected by 
over-ionization. 

Regardless of the physical reasons for the observed be- 
havior, we computed the mean [X/Fe] ratios for sodium, 
aluminum, and calcium considering only the ONC/OBlb 
stars with temperatures higher than 4500 K (Table @| . In 
the case of titanium, where some ionized lines were mea- 
sured, we list the average coming from Til and Tin. There 
is no significant difference in the average values between the 
two subgroups and the dispersion observed in the elemental 
abundances of each subgroup is smaller than the observa- 
tional uncertainties. The average [X/Fe] values are close to 
solar, with Na, Al, and Si being slightly overabundant. 

5.3. Comparison with results of other authors 

Some of our sample stars have been observed by other au- 
thors. In par ticular, one sta r, JW157=KM Ori, is in com- 
mon with the lPadgetfj (|1996f) sample. She finds for this star 
[Fe/H]=+0.14±0.18, which disagrees significantly with our 
iron abundance (A[Fe/H]=0.29 dex). Possible explanations 
include: 

i) The Padgett abundance analysis was based on fewer 

iron lines (17) than our own (40 lines; see Table HJ); 
ii) whereas we excluded very strong lines with EW > 150 
mA (see Sect. 14. 2[) . 6 among 17 lines in the Padgett list 
have EW > 150 mA (see her Table 9); 
Hi) seven lines of the Padgett list are in common with us, 
but for four of those we find different EWs of about ±15- 
20 mA. We carefully re-measured these lines confirming 
our original values. We conclude that our determination 
is likely to be more reliable. 

Both the mean [Fe/H] value for the ONC and [X/Fe] 
ratios fou nd by ourselves are in good agreement with the 
results of ISantos et al.l (|2008l ) . who derived a mean metal- 
licity of [Fe/H]=-0.13 ± 0.06 and Ni and Si abundances 
[Ni/Fe]=-0.06±0.07 and [Si/Fe]= 0.00±0.09, respectively. 
In particular, two of three stars of their sample (JW589 
and JW157) are in common with ours. There is a reason- 
able agreement between the abundances for these two stars, 
and any differences can be accounted for by the different 
line lists and cr-clipping criteria used. 

As show n in Sect. | 4. 31 the pr esent results yield a lower 
[Fc/H] than lD'Orazi et ahl(|2009t ). The main reasons for this 
discrepancy are: 

i) our new analysis is based on more suitable spectra and 
line list of Fei, which allow us to more tighly constrain 



the abundance versus \ an d versus EW trends, hence 
the T g and £ values; 

ii) our analysis of stars with T e g 4400 K is based on 
GAIA atmospheric models that are more appropriate 
for cool stars (as discussed in Appendix [Bjl ; 

Hi) critical and careful reanalysis of the CD#4 spectra ac- 
quired by D'Orazi et al. allowed us to derive the micro- 
turbulence value of the star # 673, which had been not 
inferred by these authors. 



6. Discussion 

6.1. Triggered star formation and chemical self-enrichment 

Orion is regarded as a proto-type of triggered star for- 
matio n, where star formation h as proceeded sequen- 
tially (|Preibisch fc Zinneckerl [2006') . As mentioned in the 
Introduction, this scenario predicts a peculiar chemical en- 
richment due to contamination of material ejected from 
type-II supernovae (SNIIe) originating from a first gener- 
ation of massive stars, since these are expected to contain 
the nucleosynthetic products of the stellar interior. 

In support of this view, ICunha fc Lambertl (|1992lll994ft 
found that stars in the young subgroup Id and some of the 
slightly older subgroup lc have an abundance up to about 
40% higher than the rest of their sample. They suggested 
that the enrichment resulted from the mixing of SNIIe 
ejecta f rom the lc subgrou p to the center of the Trapezium 
cluster. ISimon-DfazI ([20101 ) derived homogeneous values of 
the oxygen and silicon abundances in stars of the four sub- 
groups (OBla,b,c,d), which had a dispersion (^0.04 dex) 
smaller than the intrinsic uncertainties (~ ±0.10 dex). 

Our results indicate that low-mass stars yi eld the same 
abun dance distribution as high-mass stars (see lSim6i> Diaz 
2010). In particular, Si is the only a-element that does not 
exhibit strong evidence of being affected by NLTE on the 
basis of our data (see Fig. ITT]) and for which we obtained 
abundances for both ONC-OBld and OBlb. We find group- 
to-group dispersions of ^0.08 and ^0.01 dex in [Si/H] and 
[Si/Fe], respectively, which are smaller than our internal 
errors (of around ±0.11-0.31 and ±0.13-0.36 dex, respec- 
tively). The other elements for which we measured the 
abundance is titanium. We find for this element a differ- 
ence of 0.08 dex for <[Ti/H]> and 0.01 dex for <[Ti/Fe]>; 
this is smaller than our internal errors (around ±0.10-0.15 
and ±0.12-0.31 dex, respectively). 

We conclude that the ONC and OBlb are character- 
ized by homogeneous silicon and titanium abundances. This 
means that even if SNII explosions occurred in OBlb, at 
the OBlb-ONC distance their ejecta did not have the con- 
ditions to chemically enrich the ONC stars, dispersing the 
element over a large volume. 

As for the difference in [Fe/H] between the ONC and 
OBlb (should it be real), we note that an inhomogencity 
in metallicity (at the level of ~ 0.05 dex) within a given 
star-forming co mplex is expected in models of hierarchical 
star formation ([Emiegreenlll998f) . Owing their chaotic and 
large-scale formation process on a 1 kpc scale, the gas in 
a giant molecular cloud will have a range of metallicities 
reflecting the background Galactic gradient. We find this 
unlikely because the separation between ONC and OBlb 
(< 50 pc) is much smaller than the scale on which the 
Galactic gradient operates. 
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6.2. The metallicity of SFRs: comparison with young open 
clusters 

Our presen t ana lysis reinforces the conclusion of 
ISantos et al.l (|2008f) that none of the SFRs with available 
metallicity is more metal-rich than the Sun; the majority 
of the m are indeed slightly more metal-poor. ISantos et al.l 
(2008) suggest that, if the lower-than-solar metallicities of 
SFRs were confirmed, this would imply that either the Sun 
was formed in an inner region of the Milky Way disk, or 
that the nearby interstellar medium experienced a recent 

inf all of metal-poor gas . 

iD'Orazi &: Randichl (|2009t ) demonstrated that the abun- 
dance pattern of young open clusters in the solar neigh- 
borhood is identical to the solar distribution, concluding 
that the Sun was most likely born at the present location. 
In Fig. [T31 we compare the distribution of [Fe/H] of the 
SFRs within 500 pc of the Sun with that of i) open clusters 
younger than ^150 Myr within the same distance from 
the Sun, and ii) young nearby loose associations. In addi- 
tion to the [Fe/H] values for the ONC and OBlb derived 
here, metallicity determinations for other SFRs were re- 
trieved from Santos et al. (2008 - Chamaeleon, p Ophiucus, 
Corona Australis, Lupus), Gonzalez-Hernandez et al. (2008 
- a Orionis), and D'Orazi et al. (in preparation - Taurus). 
[Fe/H] values for the young associations were taken from 
Viana Almeida et al. (2009 - their 'uncorrected' values were 
considered) . 

The figure shows that the three distributions, in par- 
ticular that of open clusters, are characterized by a small 
dispersion; however, the cluster distribution is shifted to- 
wards somewhat higher metallicities than the SFRs and 
young associations. We obtain average values of [Fe/H] = 
-0.06 ± 0.04, -0.06 ± 0.04, and -0.01 ± 0.03 for the SFRs, 




[Fe/H] 



Fig. 13. [Fe/H] distribution for: SFRs within 500 pc 
from the Sun (dashed histogram), young nearby loose 
associations (dot-dashed line), and open clusters younger 
than 150 Myr and within 500 pc from the Sun (solid line). 
Sources of [Fe/H] val ues for the young op en clusters are the 
following: a Per s ei (iBoesgaard fc Fr iell I1990T) . NGC 2516 



(iTerndrup et al.l 120021) . NGC24 51A/B dHunsch et al 
20041) . Blancol (|Ford et al.l 120051). IC 4665 (IShen et al 
2005), IC 2602 and IC 2391 (ID'Orazi fc RandichT l2009') 



and Pleiades (jSoderblom et al.l f2009'). 



associations, and open clusters, respectively. In addition, 
none of the open clusters is as metal-poor as the most 
metal-poor SFR (the ONC) and only three out of nine SFRs 
fall within the open cluster distribution. In other words, 
not only the SFRs are more metal-poor than the Sun, but 
they are on average more metal-poor than young open clus- 
ters, which should be representative of the metallicity in 
the solar neighborhood. As for the loose associations, their 
distribution is in closer agreement with that of the SFRs 
than with the clusters; we note, however, that the 'cor- 
rected' values of IViana Almeida et al.l ([2009) would yield a 
hi gher metallici t y. Hig her values of [Fe/H] were also derived 
by iRoias et~ai1 (|2008l ) for a sample of Tucana-Horologium 
members. 

Focusing on the comparison between clusters and SFRs, 
the offset between their mean [Fe/H] values is small and 
still based on relatively small number statistics. One might 
have missed the more metal-rich SFRs or, viceversa, the 
more metal-poor open clusters in the solar vicinity. In ad- 
dition, whereas an agreement on the metallicity of the ONC 
seems now to have been reache d, discrepanc i es stil l ex- 
ist for other regions; for example, ISantos et a l. (2008) find 
[Fe/H]= — 0.08 ± 0.12 for Ophiuc us, comp a red to the value 
[Fe/H]= 0.08 ± 0.07 derived by iPadeettl (fl99l . This in- 
dicates that additional homogeneous studies must be per- 
formed, before the conclusion that SFRs are slightly more 
metal-poor than the Sun and the open clusters can be 
definitively drawn. 

With this caveat in mind, we would like to point out 
that, given the young age of the clusters, stellar migration 
is unlikely to be the reason for the difference in the [Fc/H] 
distributions of the open clusters and SFRs. The difference 
likely reflects a difference in the interstellar gas from which 
members of SFRs and young clusters formed. This in turn 
must be a relic of the process of star formation in the solar 
neighborhood, rather than an effect of chemical evolution, 
given the short timescales involved, and that in any case 
chemical evolution would lead to the younger regions (i.e., 
the SFRs) being more metal-rich than the older clusters. 

6.3. The case of P 1455: a metal-rich star? 

P1455 is significantly more metal-rich than the ONC stars. 
For this target, we derived an iro n abundance of . 1 1 ± . 09 
in very good agreement with the ICunha et al.l (|1995l ) value 
of [Fe/H]=0.08±0.15, which confirms that this star is more 
metal-rich than the other ONC targets at the 2-er level. 
Although this is still marginally consistent with what is 
expected from statistical fluctuations, some discussion of 
this object would be merited. 

At present, we do not have reasons to consider this 
star as a non-member of the ONC. Its radial velocity and 
proper motion are consistent with membership. We also 
note that we obtained an independent RV estimate using 
spectra acquired with HARPS (High Acc uracy Radial ve- 
locity Planet Searcher: iMavor et al.ll2003[) . yielding a mean 
V Iad of 21.893±0.014 km s" 1 . 

As noted in Sect. [2 this star is farther away from the 
main cluster than the other targets. Its higher metallicity 
may imply that the ONC region is inhomogeneous on scales 
of 1-2 pc and that one could search for other metal-rich stars 
close to P1455. 

We recall that whether super-solar metalli city stars ex- 
ists in SFRs has been disputed for a few years (jSantos et al.l 
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2008, and references therein), and is relevant because cir- 
cumstellar disks around young stars are the birthplace of 
planets. Thus, 'P1455-like cases' in SFRs may well be good 
targets for exoplanet searches. 

7. Conclusions 

We have presented new measurements of the abundances 
of iron-peak elements and a-elements in two subgroups 
of the Orion complex, the Orion Nebula Cluster (ONC), 
and the OBlb sub-association, derived from flames/uves 
high-resolution spectroscopy. Our main results can be sum- 
marized as follows: 

— The ONC and OBlb have mean iron abundances of 
-0.11 ± 0.08 and -0.05 ± 0.05, respectively. A likely 
member of OBla has an abundance of [Fe/H]=— 0.08 ± 
0.15. We can exclude the metallicity of Orion being 
above the solar value. 

— The ONC and OBlb are characterized by a small scat- 
ter in iron abundances, with the only exception of 
P1455 in the ONC, which we co nfirm to be metal- 
rich, as found in previous st udies (jCunha et al.lll998t 
iCavrel de Strobel et all[200l . 

— In the temperature range where NLTE effects are less 
evident (i.e. ^4500 K), there is no presence in the 
ONC sample of star-to-star inhomogeneity in the abun- 
dances of the elements strongly affected by these effects 
(namely Na, Al, Ca, Til). Owing to the lower tem- 
peratures of the OBlb sample, we are unable to draw 
any conclusions about these elements. For elements not 
strongly affected by NTLE effects (namely, Si and Ni) 
both the ONC and OBlb do not show any star-to-star 
abundance inhomogeneity. 

— The two sub-associations analyzed here have similar so- 
lar abundances of the a-elements silicon and titanium 
(the latter obtained by averaging the abundance of Ti I 
and Tin). Similar nickel abundances were found for the 
two Orion subgroups. No evidence of self-enrichment 
from OBlb to the ONC is found. 

— Star-forming regions and open clusters younger than 
150 Myr and within 500 pc of the Sun were found to 
have a small offset between their mean iron abundance, 
with the former being more metal-poor than the latter. 
Owing to the young age of both sets of stars, this offset 
probably reflects a difference in the properties of the in- 
terstellar gas from which members of SFRs and young 
clusters formed. More homogeneous studies are required 
to draw definitive conclusions. 
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Table A.l. Wavelength, elements, excitation potential, and oscillator strength of all the elements are listed. 
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Tabic FrOl (continued) 
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Tabic [A~T1 (continued) 



\ 

s\ 




X 


iu s>yj 


( A ) 




(evj 




63fo.24/ 


JNi I 


4.154 


n o on 

— O.ooO 


6384.668 


Nil 


4.154 


-1.130 


6586.308 


Nil 


1.951 


-2.810 


6598.611 


Nil 


4.236 


-0.980 


6635.137 


Nil 


4.419 


-0.830 


6767.784 


Nil 


1.830 


-2.060 


6772.321 


Nil 


3.660 


-0.960 



Appendix B: Elemental abundance analysis of cool stars: dependence on model atmospheres 

Many possible fallacies can affect the process going from the spectroscopic stellar observations to the derivation of the 
chemical composition using atomic parameters and the derivation of stellar parameters. We focus here on the role of 
model atmospheres. In particular, we show how the use of different model atmospheres leads to different results in 
metallicity and other elemental abundances (i.e., sodium, aluminum, silicon, calcium, titanium, and nickel). 

B.l. The test 

In the following, we present our starting points: 

— We considered three young members of ONC and OBlb (namely, CVS0159, CVS0118, and KM Ori) because they 
cover a wide range in effective temperature (T e g ~ 4000 — 4700 K) and surface gravity (logg ~ 3.0 — 4.5), but our 
dis cussion can be e xtend ed to all late-G/early-M stars. We also considered, as a comparison, a solar spectrum acquired 
by iRandich et al.l (|2006h with flames / uves at a similar resolution of the other spectra. 

— Abundance analysis was carried out following th e steps given in Sect. 0J 

— We considered low-resolution f 20 A) ATLAS0 (|KuruczHl993l) and high-resolution (2 A) GAIAB (|Hauschildt et al.l 
1 19991 : iBrott fc Hauschildtl 12010) synthetic spectra to evaluate the continuum flux around lines and in photometric 
bands normally used for line EW measurements (see Sect. IB.2[) . ATLAS spectra cover the ultraviolet (1000 A) to 
infrared (lOum) sp ectral range, while GAIA s pectra cover the 300 A < A < 100 /im wavelength range. 

— lKuruc3 (jl993l ) and IBrott fc Hausc hildt (2010) grids of plane parallel model atmospheres were considered for the 
abundance measurements (see Sect. IB. 3[ ). ATLAS includes atmosphere models with metallicities —5.0 <[Fe/H]< +1.0, 
gravity range 0.0 < log.g < 5.0, and 3500< T cff <10 000 K. GAIA model atmospheres span in 2000 < T off < 10000 K, 
0.0 < logg < 5.5, and —4.0 <[Fe/H]< +0.5. Model atmospheres for specific stellar parameters o f our interest wer e 
generated by interpolating in the original ATLAS and GAlA grids (see the procedure described bv lBean et al.| [2006). 

B.2. Implications on continuum flux 

One of the most important improvements made to the GAIA models was the inclusion of millions of molecular lines in the 
line list. This is of paramount importance when computing the band opacity, in addition to the line opacity. The effects 
of band opacities on the continuum flux are most pronounced in the optical domain that is largely used for abundance 
measurements (namely, 4000-8000 A). 

To help identify the range in effective temperature where the two grids of models can be used for abundance mea- 
surements, we calculated the GAIA average continuum fluxes in 20 A windows centered on A4900, 5600, 6300, and 7500 
A, which are typical regions used for abundance determinations. These fluxes were evaluated for solar-scaled chemical 
composition. Continuum flux at the same wavelengths were also considered for ATLAS low-resolution spectra (sampled at 
20 A) of solar abundance. The comparison of these fluxes is shown in Fig. IB. 11 for log g = 4.0 and for 3000 < T c $ < 7500 
K, which are a typical gravity and temperatures of low-mass members of star-forming regions. At all temperatures, 
the flux obtained with the GAIA model is lower than the flux obtained using ATLAS models, but for T cff < 4400 K 
(depending on the line wavelength) the flux decrement of GAIA spectra is more pronounced than the ATLAS spectra. 
This is particularly evident, for instance, at A = 5600, 6300 A, and is due to the formation in stellar spectra of molecular 
bands, such as metal oxide (most of all TiO, but also VO), hydroxide (such as OH), hybrids (such as CaH, FeH, MgH) 
in the visible, and CO and H2O in the near infrared. These bands are not accurately reproduced by ATLAS models, 
which lack line opacity computations for both triatomic molecules (with the exception of H2O) and numerous diatomic 
molecular transitions (such as VO). 

Since the lines used for abundance measurements are spread over wide spectral ranges (typically in the 4000 — 8000 
A range), we also calculated the synthetic fluxes in the Johnson j3]/i?J-bands by integrating the synthetic ATLAS and 
GAIA spectra, weighted by the Johnson transmission curve of the BVRI filters. The results are shown in Fig. IB. 21 for 

2 http://kurucz.harvard.edu/ 

3 http:/ /www.hs. uni-hamburg.de/EN/For/ThA/phoenix/ 
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log g = 4.0. The shift between the ATLAS and GAIA models also appears in the BVRI-Ruxes, even if it is less evident 
than the line-continuum fluxes of Fig. IB. ll because of the integration over the band wavelengths. 



•< 7 



lj- 5 



A=4900 A 



A* 



□ ATLAS 
* GAIA 



3000 4000 5000 6000 7000 

\„ (K) 



X=5600 A 



■< 7 . 

E : 



6 r 



l 5: 



□ ATLAS 



4 h 

3000 4000 5000 6000 7000 

Lf, (K) 



•< 7 

E 

O 

Cn 6 



: A=6300 A 
























X 








□ 


ATLAS : 










GAIA : 



3000 4000 



5000 6000 
\„ (K) 



7000 



•< 7 r 



Ll- 5 



4 



: \= 


7500 A 








a - B -.#-* D ** 7 




a* 




: 










□ ATLAS : 
* GAIA : 



3000 4000 



5000 6000 
T, H (K) 



7000 



Fig. B.l. Comparison between continuum flux at A4900, 5600, 6300, 7500 A obtained with ATLAS spectra (squares and 
dashed line) and GAIA spectra (asterisks and dotted line) at logg = 4.0 as a function of T c ff- The lines represent an 
interpolation through the points. 



B.3. Implications on abundance determination 

To search for the effect on abundances, we compared the mctallicitics and the Na, Al, Si, Ca, Ti, Ni abundances obtained 
using ATLAS and GAIA models for CVS0159, CVS0118, KM Ori, and the Sun. The results are listed in Table |B7H for 
both ATLAS and GAIA models. 



B.3.1. Iron 



The comparison of metallicities is shown in Fig. IB.3l where the difference between the two models increases with decreasing 
temperature, because of the presence of the mentioned bands. In particular, at T c s = 4000 — 4300 K the difference in 
abundance is ±0.07 — 0.08 dex, while for the Sun the difference is only ±0.01 dex (see also Table |B~H . As a consequence, 
for a differential abundance analysis with respect to the Sun, for KM Ori (at 4700 K) we do not find almost any 
difference between the two models, while for both CVS0159 (at 4000 K) and CVS0118 (at 4300 K) the differences are 
[Fe / H] atlas — P e /H] gaia =±0-07, ±0.08 dex, respectively. This means that the strong departure of the GAIA model 
from the ATLAS behavior at ~ 4400 K shown in Fig. IB. II affects iron abundance, leading to similar differences in stars 
with T e ff ^ 4400 K. The lower metallicity resulting from GAIA models (with respect to the ATLAS models) depends on 
the more signidicant formation of molecules in atmospheres with lower temperatures. The molecular opacity considered 
in the GAIA models indeed leads to a redistribution of the flux, which is on average lower than the ATLAS one, because 
of the molecular absorption (line blanketing). Lower flux yields lower intrinsic line equivalent widths, which can be 
reproduced by lower iron abundances. 
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Fig. B.2. Comparison between continuum flux at the Johnson f?Fi?/-bands obtained with ATLAS spectra (squares and 
dashed line) and GAIA spectra (asterisks and dotted line) at log g — 4.0 as a function of T c g. The lines represent an 
interpolation through the points. 



B.3.2. Other elements 

In Table IB. It we summarize how the use of different models can affect the elemental abundance. Here, the comparison 
between ATLAS and GAIA grids is listed for Na, Al, Si, Ca, Ti, and Ni. While lines of elements across over the whole 
spectrum (~ 4800 — 6800 A) infer very different results (such as Ni, besides Fe), elements such as Al, with only two lines 
at 6696 A and 6698 A close each other and not strongly affected by band opacity, lead to similar GAIA and ATLAS 
abundances. 



B.4. Concluding... 

We find T e g « 4400 K to be the lower limit where the models in which the line opacity computations are not fully 
treated, such as ATLAS, can be applied in an abundance analysis. This has been demonstrated for both iron abundance, 
typically derived by many lines, and other elements (a- and iron-peak elements) typically used as tracers of chemical 
enrichment. 
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Fig. B.3. Comparison of iron abundances derived by using ATLAS and GAIA model atmospheres as a function of 
equivalent width (EW) and line excitation potential (EP). These examples display the results obtained for stars with 
four different temperatures. Circles and asterisks refer to abundances derived with ATLAS and GAIA models, respectively, 
while solid and dashed lines represent their mean values. Vertical thick and thin bars in the logn(Fe) vs. EW panels are 
the standard deviations around the average iron abundances. 
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Table B.l. Examples of Fe, Na, Al, Si, Ca, Ti, and Ni mean abundances obtained for stars in the Orion complex using 
ATLAS (A) and GAIA (G) models. 



Name 


CVS0159 


CVSOH8 


KM Ori 


Sun 


T e ff 




4000 K 


4300 K 


4700 K 


5770 K 






O O 

o.y 


A O 


O 1 

o.l 


A A A 

4.44 


< logn(Fe) A 


> 


7.46 ± 0.24 


7.60 ±0.11 


7.37 ±0.07 


7.52 ±0.02 


< logn(Fe)G 


> 


*7 on _l n oo 
/.39 ± 0.26 


*7 co _i_ n i o 
1 .01 ± 0.12 


*7 O A _1_ H AT 

(.34 ± 0.0/ 


^ n i n no 
l.oi ± 0.U2 


< logn(Na) A 


> 


6.04 ±0.07 


6.11 ±0.03 


6.28 ±0.05 


6.31 ±0.04 


< logn(Na)G 


> 


a no _l n r\a 
6. 08 ± 0.06 


a i i _i_ n no 
6.11 ± 0.03 


a _i_ n nt; 
6.26 ± 0.U6 


p. on _i_ n no 
6.29 ± 0.U3 


< logn(Al) A 


> 


6.27 ±0.04 


6.41 ±0.02 


6.45 ±0.12 


6.48 ±0.03 


< logn(Al) G 


> 


6.26 ± 0.03 


g oo _j_ n m 
6.38 ± 0.01 


a ao _i_ n i a 
6.4z ± 0.14 


p. a*7 _i_ n no 
6.4/ ± 0.U3 


< logn(Si)A 


> 




7.65 ±0.09 


7.55 ±0.10 


7.56 ± 0.03 


< logn(Si)G 


> 




7.65 ±0.09 


7.51 ±0.11 


7.53 ± 0.03 


< logn(Ca)A 


> 


6.10 ±0.04 


6.17 ±0.11 


6.22 ±0.10 


6.35 ±0.03 


< logn(Ca)G 


> 


6.12 ±0.03 


6.16 ±0.11 


6.17±0.10 


6.34 ±0.03 


< logn(Tii) A 


> 


4.45 ± 0.23 


4.70 ±0.07 


4.77 ±0.08 


4.97 ±0.02 


< logn(Til) G 


> 


4.49 ± 0.22 


4.68 ± 0.08 


4.74 ± 0.09 


4.97 ±0.02 


< logn(Ni) A 


> 


6.20 ±0.03 


6.30 ±0.11 


6.01 ±0.10 


6.26 ±0.03 


< logn(Ni) G 


> 


6.11 ±0.03 


6.22 ±0.10 


6.02 ±0.11 


6.24 ±0.02 
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